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Abstract 

Motivated by experimental studies on the anomalous diffusion of bi- 
ological populations, we introduce a nonlocal differential operator which 
can be interpreted as the spectral square root of the Laplacian in bounded 
domains with Neumann homogeneous boundary conditions. Moreover, 
we study related linear and nonlinear problems exploiting a local realiza- 
tion of such operator as performed in 1 7 for Dirichlet homogeneous data. 
In particular we tackle a class of nonautonomous nonlinearities of logistic 
type, proving some existence and uniqueness results for positive solutions 
by means of variational methods and bifurcation theory. 

1 Introduction 

Nonlocal operators, and notably fractional ones, are a classical topic in har- 
monic analysis and operator theory, and they are recently becoming impres- 
sively popular because of their connection with many real-world phenomena, 
from physics |2Ul[T4j|2I] to mathematical nonlinear analysis 1 1 , 24 1 , from finance 
IH [T3] to ecology (6l [23j [T71 15] . A typical example in this context is provided by 
Levy flights in ecology: optimal search theory predicts that predators should 
adopt search strategies based on long jumps -frequently called Levy flights- 
where prey is sparse and distributed unpredictably, Brownian motion being more 
efficient only for locating abundant prey (see |25l [TT)). As the dynamic of 
a population dispersing via random walk is well described by a local operator 
-typically the Laplacian- Levy diffusion processes are generated by fractional 
powers of the Laplacian (-A) s for s e (0, 1) in all U N . These operators in IR^ can 
be defined equivalently in different ways, all of them enlightening their nonlocal 
nature, but, as shown in |8| and |9|, they admit also local realizations: the frac- 
tional Laplacian of a given function u corresponds to the Dirichlet to Neumann 
map of a suitable extension of u to U N x [0, +00). On the contrary, on bounded 
domains, different not equivalent definitions are available (see e.g. H5l|3J[7| and 
references therein). This variety reflects the different ways in which the bound- 
ary conditions can be understood in the definition of the nonlocal operator. In 



particular, we wish to mention the recent paper by Cabre and Tan 1 7 1 , where the 
operator (-A) 1/2 on a bounded domain Q e U N and associated to homogenous 
Dirichlet boundary conditions is defined by Fourier series, using a basis of corre- 
sponding eigenfunctions of - A. Their point of view allows to recover also in the 
case of a bounded domain the aforementioned local realization: indeed, inter- 
preting D. = Q x {0} as a part of the boundary of the cylinder O x (0, +oo) a U N+l , 
the Dirichlet spectral square root of the Laplacian coincides with the Dirichlet to 
Neumann map for functions which are harmonic in the cylinder and zero on its 
lateral surface. These arguments can be extended also to different powers of - A, 
see Q2] . On the other hand, in population dynamic, Neumann boundary data 
are as natural as Dirichlet ones, as they represent a boundary acting as a perfect 
barrier for the population. The aim of this paper is then to provide a first contri- 
bution in the study of the spectral square root of the Laplacian with Neumann 
boundary conditions. 

Inspired by Q, our first goal is to provide a formulation of the problem 



(1.1) 



|(-A) 1/2 a = / inQ, 
\d v u = on3Q, 



where O is a C 2,a bounded domain in R , N > 1, and / can be thought, for 
instance, as an L 2 (Q) function. To this aim, let us denote with {</>fc} fc>0 an or- 
thonormal basis in L 2 (Q) formed by eigenfunctions associated to eigenvalues 
file of the Laplace operator subjected to homogenous Neumann boundary con- 
ditions, that is 



I -A0fc = jUjfc0jfc inQ, 
\d v (p k = on3Q. 



(1.2) 

We can define the operator (-A) 1/2 : H 1 (Q) — I 2 (Q)) by 



(1.3) (-A) u= YlL 1 ic u k<Pk forugivenby u=^ui c (f>j c . 

k=l k=0 



The first series in (L3) starts from k = 1 since the first eigenvalue and the cor- 
responding eigenfunction in (1.2) are given by (jUo></>o) = (0, 1/ vTOl)- This sim- 
ple difference with the Laplacian subjected to homogeneous Dirichlet boundary 
conditions has considerable effects. First of all, this implies that (-A) 1 ' 2 , as the 
usual Neumann Laplacian, has a nontrivial kernel made of the constant func- 
tions, then it is not an invertible operator and cannot be solved without 
imposing additional conditions on the datum /; on the other hand, given any u 
defined on Q, its harmonic extension on f := fl x (0, +oo) having zero normal 
derivative on the lateral surface needs not to belong to any Sobolev space, as 
constant functions show. These features has to be taken into account when es- 
tablishing the functional framework where to set the variational formulation of 
{1.1) . In this direction, we will first provide a proper interpretation of (LTJ , and 
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a corresponding local realization, in the zero mean setting. To this aim, let us 
introduce the space of functions defined in the cylinder ^ 

je x { c €'\:=\v^B x { c €'): f v{x,y)dx = 0,Vye[0,+oo)\. 
An easy application of the Poincare-Wirtinger inequality shows that we can choose 



2.2 



and 



as a norm of v e J€ 1 (Sf ) the L 2 norm of the gradient of v (see Proposition ; 
Lemma |23) . It comes out that, when the datum / has zero mean, a possible so- 
lution of (1.1) is the trace of a function belonging to J€ 1 (^). The corresponding 
space of traces can be equivalently defined in different ways, since Proposition 
IO] shows that 

^ lt2 (D.):=\ueH lt2 (D.): J u{x) dx = oj = {u = v{x,0) : ye^f 1 ^)} 

{+oo +oo I 

uel 2 (fi): u= Uk<pk such that £ /ij. /2 u| < +00 > . 
k=i k=i J 

In proving this result, one obtains that every u e J€ l (p.) has an harmonic ex- 
tension v e 36 1 (S^) given by 



(1.4) v[x,y) = Y, u k (p k {x)e ^* y , for (/u fc ,0 fc ) solving (l~2"l , 

fc=i 

and which is also the unique weak solution of the problem 

Ay = in^, 

(1.5) {d v v = onffix(0,+oo), 

v[x,0) = u[x) on D.. 

Thus, given u £ we can find a unique f e J€ 1 (^) solving (T3J , for which 

it is well defined the functional acting on H ll2 {D.) as 



<■ 



-d y v{-,0),g) := J Vi/-Vgdxdy, 



where g is any H 1 ^) extension of g. Since this functional is actually an ele- 
ment of the dual of J? ll2 {n), it is well defined the operator L\/2{u) = -d y v{-,0) 
between J€ 1/2 (fi) and its dual. Thus, restricting the study to the zero mean func- 
tion spaces, and taking into account equations (L3) and (1.4) , we have that Z4/2 
conincides with (-A) 1/2 , but it is invertible: for every / in the dual space of 
.^f 1/2 (Q) there exists a unique u e ^f 1/2 (Q) such that L1/2W = /, and this func- 
tion u is the trace on Q of the unique solution v e J€ Y of the problem (see 
Lemma |2T4} 

(Ay = inOxfO.+co), 
d v v = on50x(0,+oo), 
d v v(x,0) = f(x) on D.. 
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The link between Z4/2 and (- A) ' now becomes transparent since 



(-A) 1/2 (w) = L 1/2 



that is, the image of a function u trough (-A) 1/2 is the same of the one yield by 



L1/2 acting on the zero mean component of u (see Definition 2.12 1. In this way 
we have recovered the local realization of (-A) 1/2 as a map Dirichlet-Neumann 
since 

u - + u\ = -d y v{x, 0) = d v v 



(-A) 1/2 u = L 1/2 



where v solves (1.5) with Dirichlet datum u = u- f u instead of u. Therefore, if 
/ has zero mean, denoting with u{x) = v[x,0) the unique solution of Z4/2W = / 
then the solutions set of is given by u + h for heU. 

Since we are interested in ecological applications, as a first study we focus our 
attention on the logistic equation. More precisely, consider a population dis- 
persing via the above defined anomalous diffusion in a bounded region £1, with 
Neumann boundary conditions, growing logistically within the region; then u, 
the population density, solves the diffusive equation 

u t = d{-\) ll2 u{x, t) + u(x, t)(m(x) - u(x, t)) in O. x (0, +00), 
d v u = on d£l x (0, +00), 

u{x,0) = uq{x) inO, 

where d > acts as a diffusion coefficient, the term - u 2 express the self-limi- 
tation of the population and m e C°' l [D.) corresponds to the birth rate of the 
population if self-limitation is ignored. The weight m may be positive or nega- 
tive in different regions, denoting favorable or hostile habitat, respectively. The 
stationary states of this equation are the solutions of the following nonlinear 
problem 



(1.7) 



[ (-A) 1/2 k = Aw(m(x) - u) inf2, 
\d v u = Q on dQ., 



where A = lid > 0. When the diffusion follows the rules of the Brownian motion 
this model has been introduced in |26| and studied by many authors (see HU 
and the references therein). One of the major task in this problem is describing 
how favorable and unfavorable habitats, represented by the interaction between 
A and m, affects the overall suitability of an environment for a given populations 
flUl . The typical known facts for the stationary problem associated to Brownian 
motion can be summarized as follows: 

Theorem 1.1 (|16, 28|). i) If the function m has negative mean inside Q and it 
is positive somewhere, then there exists a positive number Ji l such that for every 
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A > p 1 there exists a unique positive solution u\ 

{- Aw = Xu(m(x) - u) inn, 
o v u = onoD., 

and ux^ as A — >• ^ . 

ii) Ifm has nonnegative average, then for every A > there exists a unique positive 
solution U\ of (L8) and u\ — h* as A — + , /or h* expressed by 

(1.9) h* = £ m{x)dx - — [ m[x)dx. 

J I d| Jn 

The number ju 1 appearing in i) is the first positive eigenvalue with positive eigen- 
function of the operator - A with Neumann boundary condition and with a weight 
m satisfying the hypotheses in i) . 

In our situation, we have, first of all, to clarify that by a weak positive solution 
of (L7) we mean a function u e H 1 (f2), u{x) > 0, u{x) = ii(x) + h with h eU + 
and u e ^ 1/2 (Q), so that u{x) = v{x,0) for V e ^f 1 ^) and (V, h) e Je l &) x U is 
a weak solution of the nonlinear problem 



(1.10) 



A£ = inQx(0,+oo), 
d v i> = on x (0, +oo), 

d v v = \{v + h)[m{x) - v - h) onQx{0}, 

/ X{v{x, 0) + h) {m{x) - v{x, 0) -h)dx = 0, 
Jn 



in 

in the sense that 

/ VvVy/dxdy=l X{v + h){m[x)-v-h)y/dx Vy/ e ^f 1 ^), 
J^ Jn 

\ X(_v + h)(m(x) - v - h)dx = 0. 
Jn 

In other words, we impose that the right hand side has zero mean, choosing, in 
this way, the mean of a solution u as h = h u . Then we obtain the well posedeness 
of the problem 

{L\i2U = X{ii + h u ){m(x) - u-h u ) infi, 
d v u = ondQ, 

since now the right hand side has zero mean, and we obtain in this way the zero 
part mean of u. Moreover, notice that the mean of the function v{x,y) solution 
of (LIP) with v (x, y) = v + h u and v(x, 0) = is exactly the mean of u. 
Our main existence result is the following 

Theorem 1.2. Let m e C 0,1 (Q), m^0. Then the following conclusion hold: 

i) If f n m{x)dx < and there exists Xq e Q such that m(xo) > 0, then there exists 
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a positive number A\ such that for every A > Ay there exists a unique positive 
solution ux of (L7) , with u\ —- as A — ► A* . 

ii) If f n m[x) dx > 0, then for every A > there exists a unique positive solution ux 
of (L7) with u\ -> ft* /or A — + . 

As in the standard diffusion case, A\ is the first positive eigenvalue with positive 
eigenfunction of the problem 

J(-A) 1/2 w = Am{x)u inQ, 
|3 v u = on3Q. 



which existence is proved in Theorem 3.7 Theorem 1.2 will be obtained via clas- 
sical bifurcation theory, indeed, in case i), we can show that a smooth carte- 
sian branch of positive solutions bifurcates from the trivial solution (A, h, u) = 
(Ai,0,0), this branch can be continued in all the interval (Ai,+oo), and contains 
all the positive solutions of {1.7) , that is to say that for every A > there exists 



a unique positive solution (see Proposition 3.15 and Theorem 3.20 . We tackle 
case ii) first assuming that the mean of m is positive. This allows us choose as a 
bifurcation parameter h, the future mean of u, instead of A, and find a branch bi- 
furcating from the trivial solution (A, h, u) = (0, h*,0), with h* defined as in (L9) . 
As in the previous case we can show that this branch is global and contains all 



the positive solutions (see Proposition 3.16 and Theorem 3.20}. Finally, we com 



plete the proof of case ii) by approximation in Theorem 3.22 



All the effort made in finding the proper formulation for the linear and the non- 
linear problem enables us to prove the existence results for {1.7) , which are in 
accordance with the case of standard diffusion. But, trying to enlighten the dif- 
ferences between the two models, one has to take care of the eigenvalues ap- 
pearing in Theorems j I • I l anc ^ ^ m &t is Ji^ and Ai. Since such eigenvalues act as 
a survival threshold in hostile habitat, it is a natural question to wonder which is 
the lowest one, indeed this indicates whether or not the fractional search strat- 
egy is preferable with respect to the brownian one. This appears to be a difficult 
question, since the eigenvalues depend in a nontrivial way on m, and also on the 
sequence {pk)k defined in {1.2) . At the end of Section[3]we report some simple 
numerical experiments to hint such complexity. 



2 Functional setting 

In this section we will introduce the functional spaces where the spectral Lapla- 
cian associated to homogeneous Neumann boundary conditions will be defined. 
Moreover, we will study the main properties of this operator and find the proper 
conditions under which the inverse operator is well defined. Finally, we will 
prove summability and regularity properties enjoyed by the solutions of the lin- 
ear problem. 

Throughout the paper Q is a C 2,a bounded domain and we will use the notation 
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^ = (lx (0,+oo). 

In this plan we will make use of the following projections operators. 
Definition 2.1. Let us define the operators : H 1 {<€) — H l {<€) by 

Acgv:= I v(x,-)dx= — [ v{x,-)dx, Zgv:=v-A<gv, 
Jn Jn 

for \Q\ denoting the Lebesgue measure of the domain Q. A^ and Z<g> give the 
average (with respect to x) and the zero-averaged part of a function v, respec- 
tively. Analogously, for u e H 1I2 {Q), we write 

(2.1) Aau:= f u(_x)dx, ZqW := u - Anu. 

Jn 

When no confusion is possible, we drop the subscript in A, Z. 

It is standard to prove that, in both cases, A and Z are linear and continuous, 
and that trn o Z% = Zn° tr^. Since the integration in the definition of Ag> is per- 
formed only with respect to the x variable, it is natural to interpret the image of a 
function v through the operator A^ as a function of one variable. A^ v{y) enjoys 
the following properties. 

Proposition 2.2. If v e H 1 ^) then A^v e ff^Oj+oo). In particular, it is a con- 
tinuous function up to + , and it vanishes as y tends to infinity. 

Proof. Since d y v[-,y) e L 2 (O) for almost every y, we can compute (A<# v)'{y) and 
obtain, by Holder's inequality, 

n+oo r+oo i(r \ 2 

Jo [iAegvy) dy= L wAL dyVdx ) dy 

*rMa [dyv{2dx ) dy<+o °- 

As a consequence, A<g v e H l (0, oo) , so that it is continuous in y and it vanishes 
as y tends to +oo. □ 

Introducing the following functional spaces 

je l ffl:=YszA< € = \vs.H x {<€)•. j v{x,y) dx = 0, Vy e (0, +oo)| , 



(2.2) 



je ll2 (n):=KerA a = \ue H ll2 (D.): \ u{x)dx = 0\, 

Jn 



it is worth noticing that the former is well defined by Proposition 2.2 Moreover, 
we can choose as a norm on 36 1 {*€) the quantity 

(2.3) lli/||^ w :=l|V< 2 

as it is equivalent to the H 1 -norm thanks to the following lemma. 
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Lemma 2.3. There exists a positive constant K such that for every v e J€ l i^€) it 
holds 

\\v\\ L 2<K\\Vv\\ L 2. 

Proof. We set V x v = [d Xl v,..., d Xn v) and we notice that for any v e J€ Y ["io) the 
Poincare-Wirtinger inequality implies 

n+oo n n+oo I n \ 

\\v\\ 2 L 2 = J o dy J^v 2 {x,y)dx< \c pw \V x v{x,y)\ 2 dx\ dy 

^c pw \ \V x v{x,y)\ 2 dxdy<c pw l \Mv{x,y)\ 2 dxdy = K 2 \\Slv\\ 2 2 

proving the claim. □ 

The following proposition gives a complete description of the space J€ Y,1 {V£). 

Proposition 2.4. Lef^f 1/2 (Q) he defined in l |2.2| l. Then the following conclusions 
hold: 

(i) je ll2 {D.) = {u = tr a v: ye^f 1 ^)} 

{+00 +00 I 

u e L 2 (Q) : u = u^ip^ such that £ p^u 2 ^ < + °° \ > 
k=l k=l J 

(ii) je ll2 {Zl) is an Hilbert space with the norm 

\ r r \u{x)-u{x')\ 2 , , ,\ m 

m ^ = \LL \x- x >\^ dxdx \ 

equivalent to the usual one in H ll2 {Q). 

Proof. Since Q is of class C 2 ' a , we have that H 1I2 (D.) can be equivalently charac- 
terized as {u = tinv : v e H 1 ^)}, where we write ixv = v\n = v{-,0). Then Propo- 
sition [Z2| provides the inclusion 

{u = tr n v. ye^mfc^p. 

In order to show the opposite one, consider u e J€ (Q) and consider ve H 1 
such that u = tr^ v. Notice that v e JC 1 and Proposition |2.2| implies that 



tx a {Z^v) = ti a (v-A<gv) = u-f v{x,0)dx= u 

Jn 

then we have found v = Z^v belonging to je 1 ^) and such that u = tin(v), 
yielding the first equality in (i). As far as the second equality is concerned, we 
start by proving the inclusion 

{+00 +00 I 

ue L 2 (Q) : u= £ U]c(p]c such that £ M^ 2 " 2 ; < +°° \ ■ 
k=l fc=l J 
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Indeed any v e 36 1 { C €) can be written as v(x,y) = Lfc>i ffc(y)</>fc(x), with 

dy 



r+oo I 

MjPM= I \Y.VkV k {y) 2 + v' k {y) 
then 



X Mfc^fcfy) 2 < +°° a.e. in (0, +00). 

fc>i 

Let us fix y such that £ Mjt'^fcCy) 2 i s finite and take u = tr^v = £ ffc(O)0fc(Jc). 

fc>i fc>i 

We have 



2 



£ ( y 2|^ /2 ^(y)^(y)|rfy; 



fc>i fc>i 



implying the desired inclusion. On the other hand, let Xfc>i M^' 2 M f < + °°» ano " 
let us define 

+ °° 1/2 

(2.4) y(Jc,y)= E u k <p k {x)e-^ y . 

k=l 

It is a direct check to verify that 1/ £ <y5f 1 (9?) (see also Lemma 2. 10 in 1 7 1) , obtain- 
ing that all the equalities in (i) hold. 

Let us now show conclusion (ii), starting with proving that there exist constants 
A, B such that 

(2.5) A||«|| H i/ 2(n) < || u\\jfii2 m < B\\u\\ H v2 {n] 
As 

2 2 r r |u(x)-u(x')l 2 . , , 

l«l flinw) = l»«l tf+ JJ o |x _ x/|W+1 dxdx. 

The right hand side inequality holds for B = 1; in order to show the left hand 
side inequality, let us argue by contradiction and suppose that there exists a se- 
quence u n £ 36 1I2 {D.), with ||u ra |li,2(Q) = 1 and ||u w ||^i/2(0) — 0. Then u n is uni- 
formly bounded in H 1/2 (Q) and there exists u e 36 112 [D.) such that u n converges 
to u weakly in H 1I2 [D.) and strongly in L 2 [D.) (notice that we do not know that 
the quantity || • \\j%>m(Q) is a norm on 36 ll2 {£L)). As a consequence, || ^ II z, 2 (O] = 1 
and 

\u{x) - u{x')\ 2 



// 

JnJn 



InJn \x-x'\ N+1 



-dxdx' = 0, 



which is an obvious contradiction. As a byproduct of inequalities ((23} we obtain 
that || • ||^i/2(Q) is a well defined norm and since 36 ll2 {£l) is a closed subspace of 
H 112 [D.) with respect to the usual norm conclusion (ii) holds. □ 

Carefully reading the proof of the second equality in (i) of the previous proposi- 
tion, one realizes that for any u e 36 II2 {$1) we can construct a suitable extension 
v e 36 1 which is harmonic and that can be written in terms of a Fourier ex- 
pansion as shown in (2.4| l . In the next lemma we provide a variational character- 
ization of such extension. 
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Lemma 2.5. For every u e JC 112 (O) there exists an unique v e J€ l (<^) achieving 

inf|||y||^ ? i (V) = J \Vv(x,y)\ 2 dxdy: v eje 1 ^), v{-,0) = u|. 
Moreover, the function v is the unique (weak) solution of the problem 

{Ay = in^, 
d v v = on dO. x (0, +00), 

v(x,0) = u{x) onD.. 

Finally 

+00 +00 

(2.7) if m(x) = ^ Ufctpidx) then w(i,y)= ^ u^Me"' 1 ' ' 

fc=i fc=i 

Proof. We observe that the functional to be minimized is simply the square of 
the norm in and the set on which we minimize is non empty and weakly 

closed thanks to the compact embedding of H ll2 {D.) in L P (Q), for any exponent 
p < 2N/{N- 1). The strict convexity of the functional implies the existence and 
uniqueness of the minimum point. 

As usual, the unique minimum point v satisfies the boundary condition on £1 
(in the H 1 ' 2 -sense) by constraint, and 

/ Vv-Vy/ = Vi/ze^? 1 ^) s.t. y/(x,0) = 0. 

As a consequence, for every ( e H 1 ^) such that ([x,0) = 0, it is possible to 
choose y/:=( - as a test function in the previous equation. This provides 

= f Vv-V(- f d y v{x,y){A^Q\y)dxdy 

n r +00 

(2.8) = / Vy-Vf-mi / {A^v)'{y){A^()'{y)dy 

Jo 

= f Vv-V( VCeH 1 ^) s.t. C(x,0)s0. 

In a standard way this implies both that v is harmonic in ^ and that it satisfies 
the boundary condition on 3Q x (0, +00) (in the if~ 1/2 -sense). 
Finally, if u{x) is given as in \2.7\ , then v as in \2.7) solves problem JZ6) and the 
uniqueness of the solution provides the claim. □ 



Definition 2.6. We will refer to the unique v solving fl2.6| l as the Neumann har- 
monic extension of the function u. 

Remark 2.7. As we already noticed, 

H V2 {n) = I u e I 2 (n) : || w|| 2 2 + [ [ |M(x) ~ M ( x)l dxdx ' < +00 1 
[ JnJn |x-x'r +i ) 

= {u = tr n v: veH 1 ^)}. 
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Furthermore, it is well known that the two norms 



JnJn \x-x'\ N ^ 



|w||^i/2 (n)|2 = inf{lkl| 2 7 i ( ^ ) : veH 1 ^), v(-,0) = wj 



are equivalent. Reasoning as in the proof of Proposition |2.4[ and taking into ac- 
we obtain that j¥ 1I2 {Q.) can be equipped with the equivalent 



2.5 



count Lemma 
norms 



2 _ r r \ujx)-u{x n 
m ^ v w JJ Q \ x - x >\n + i axax ' 



it=i 

where the terms u/ c are the Fourier coefficients of u. In particular, the harmonic 
extension of u depends on u in a linear and continuous way. 

In order to introduce and study the dual space of JF 1I2 {Q.) let us first introduce 
the following space. 

Definition 2.8. Let us define the following subspace of H~ 112 {Q). 

^f" 1/2 (D) := {/ e fT 1/2 (Q) : </, 1) = 0}, 
where <•, •> denotes the duality pairing. 

The subspace just introduced as a strict connection with the dual space of J€ 1/2 (Q) 
as well explained in the following proposition. 

Proposition 2.9. It holds je l 12 (Vt)* = je~ ll2 (P). 

Proof. We can exploit the splitting H V2 {D.) = ^f 1/2 (Q) © U in order to obtain 

je ll2 £l)* =H- V2 {n)/ ~ where /i~/ 2 ^ fiUv Hn) = f 2 \jem m . 

More precisely, on one hand if g e 1/2 (O) * then, for every c E R, 

/:=goZ a + cA n eH- 1/2 (n); 

on the other hand, if / e 1/2 (Q) then g := /|^i/2 (n) e ^ 1/2 (Q)* and 

/ = goZ n + </,l>A Q . 

Moreover, both the maps defined above are linear and continuous. This proves 
that ^? 1/2 (Q)* is isomorphic to {/ e £T 1/2 (Q) : </, 1> = c}, for every fixed c, and 
in particular for c = 0. □ 
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As a first step to arrive to a correct definition of the half Laplacian operator, let 
us prove the following lemma 

Lemma2.10. LetuE JC 1I2 {£1), andletve J€ l &} denote its Neumann harmonic 

extension. Then the functional d v v\ = -d v v(-,0): H ll2 (Q.) —> R is well de- 

lnx{0} y 

fined as 

{-d y v{-,0), g) := I Vv-Vg dxdy, 

where ge H 112 (Q) and g is any H 1 -extension ofg. Moreover, 

-d y i;(-,0)e^r 1/2 (fi). 

Proof. The functional is well defined, indeed if gi and gz are two extensions of 
g we have that (g2 - gi)ix, 0) = and, arguing as in equation < |2.8| , yields 

/ Vv-Vg2dxdy- I Vv-Vgidxdy= I Vv ■ V{g2~ gi) dxdy = 0. 

Moreover -d v v[x, 0) is linear and continuous: indeed, let us choose as an exten- 



sion of g G := Ang + g, where g is the harmonic extension of Zng; by Remark 2.7 
applied to g we have that 



\(-d y v(;0),g)\ 2 < \\v\\% lm (\A n g\ 2 + |[g 



^l|2 



<c(||g|| 2 2 + ||g||^ 1/2(n) ) = C||g| 



2 

H V2 (Q.)- 



As a consequence —d y u{x,0) e H 1/2 (Q). Finally since w{x,y) := (l-y) + belongs 
to H 1 ^€), by definition we obtain that 

(-d y v(-,0),l) = / Vv-Vwdxdy = - I v y dxdy- \ [-v{x,l) + v{x,0)]dx, 

' J<<fn{y<l} Jn 

which vanishes because v e J€ l □ 

Remark 2.11. If the harmonic extension v is more regular (for instance H 2 
then we can employ integration by parts in order to prove that the definition of 
-dyV{x,Q) given above agrees with the usual one. 

Thanks to the previous lemmas, we are now in a position to define the fractional 
operators we work with. 

Definition 2.12. We define the operator L V2 : Je ll2 {p.) — ^? _1/2 (Q) as 
(2.9) L V2 u= -dyV{-,Q). 



where v is the harmonic extension of u according to \2.6\ . Analogously, we de- 
fine the operator (-A j/) 112 : H ll2 {D.) -+ £T 1/2 (Q) by 

(-A^) 1/2 = L 1/2 °Z n . 
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In Definition 2.12 we have introduced the fractional Laplace operator associa- 
ted to homogeneous Neumann boundary conditions as a Dirichlet to Neumann 
map. Moreover, thanks to the equivalences of Proposition 2.4 we realize the 
spectral expression of this operator as explained in the following remark. 

Remark 2.13. Sin ce th e harmonic extension operator u >—■ v is linear and con- 
tinuous by Remark 2.7 we have that both Z4/2 and (- A ,jy) 112 are linear and con- 

+00 

tinuous. Moreover, if u e H^iQ) and u[x) = £ u^cp/dx), we can use equation 

fc=i 

\2.7\ to infer that d v v (x, 0) e L 2 (Q) . This allows to write 



{-Ajs) 1I2 u{x) = Li/ 2 (w)(x) = -d y v{x,0) = £ n]^ 2 u k (p k {x). 

fc=i 

In particular, ifueH 2 (fi) then 

(-A^) 1/2 o(-A^) 1/2 w = -A^w 

provides the usual Laplace operator associated to homogeneous Neumann boun- 
dary conditions on dO. 

We remark that we can think to Z4/2 as acting between J6' ll2 {D.) and its dual 
thank to Propositionlzi) While (-A j/) 112 : H 1I2 {D.) — £T 1/2 (Q) is neither injec- 



tive nor surjective, we have that Z4/2 : J6' ll2 {D.) -» J€ 1/2 (Q) is invertible. 
Lemma2.14. For every f e t ^f" 1/2 (Q) there exists a unique v z. Jt? 1 suchthat 

(2.10) f Vi/(jc,y)-V^(JC,y)djcrfy = </,^0,0)> V^e^ 1 ^). 

Moreover, the function v is the unique (weak) solution of the problem 



(2.11) 



Ay = in^, 

d v v = on dQ. x (0, +00), 

d v f (x,0) = /(x) onD.. 



Proof. The existence and uniqueness of v follows from Riesz representation The- 
orem. The fact that v satisfies l |2.11| l follows once one shows that equation (2.10) 
holds also for every y/ e H 1 ^). This can be readily done exactly as in the proof 
of Lemma 2.5 Further, it is also a consequence of the next result. □ 

The choice of JC 1 ^) as test function space is not restrictive as the following 
lemma shows. 



Lemma 2.15. Let f e J€ 1/2 (Q), and v e J? 1 ^) be defined as in Lemma 
Then: 



2.14 
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(i) there exist positive constants C, k depending on f such that, for every y > 0, 

[ \Vv{x,y)\ 2 dx<Ce~ ky ; 
Jn 

(ii) equation ^2.10) holds for any yr e Hi &~) admitting a constant C' such 
that, for every y, 

(2.12) \m;y)\\ma)^C. 

Proof. Let yr e H l (O x (0, y)). Since D. x (0, y) is bounded, we can test (27TTJ with 
yr and use integration by parts in order to obtain that, for a.e. y, 

(2.13) I Vv[x, t)- Vyr[x, t)dxdt= (f,yr[-,0)) + \ v y {x,y)yr{x,y) dx. 

J^n{t<y} Jn 

As far as the first statement is concerned, the above equation used with yr = v 
gives 

r+oo r 

$(y):= / I \Vv[x, t)\ 2 dxdt 
Jy Jn 

- \ \Vv{x,t)\ 2 dxdt- I \Vv{x,t)\ 2 dxdt 

JV J^n{t<y} 

= </,!/(-, 0)> -<j,v{; 0)>- [ v y ix,y)v{x,y)dx. 

Jn 

Then <5 is absolutely continuous and 

O(y) = - J v(x,y)Vy(x,y)dx<lj v 2 {x,y)dx^j ij v 2 (x,y)dx^ 



(n \ 1/2 / r. \ 1/2 

JjV x v(x,y)\ 2 dx\ I J v 2 {x,y)dx 

^ [ \Vv(x,y)\ 2 dx = -k(f>'{y), 
2 Jn 



c 



which implies <I>(y) < <E>(0) • e~ ky and the required inequality. 
Now we turn to the second statement. If y/ is as in its assumption then l |2.13) 
holds. In order to conclude we must prove that the last term in that equation 
vanishes as y -» +oo. But this is easily proved by applying Holder inequality and 
using the first part of the lemma. □ 

Remark2.16. In particular, the previous proposition implies that equation ^2.W\ 
holds for any yr e H l On the other hand, from its proof one can deduce that 
more general test functions are admissible, for instance functions such that their 
L 2 [D.) norm does not grow too much with respect to y. 

We are now in the position to define the inverse operator of Z4/2. 
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Definition 2.17. We define the operator T V2 : 36~ U2 {n) — 36 U2 {D.) by 
(2.14) r 1/2 (/) = frQ(i;) = i/(x,0) 

where v solves \2.1l\ . 

We collect in the following proposition the properties of 7i/ 2 . 

Proposition 2.18. The operator T\ 12 defined in \2.\A) is linear and such that L\ 12° 

Tl/2 = 7l/2°Ll/2 = Id. 

Moreover T 1/2 : £6 2 (O) := {/ £ I 2 (O), : / Q /(x) dx = 0} — «Sf 2 (Q) is compact, posi- 
tive, self-adjoint and 7\/2° Z1/2 = (—A ./k) -1 !^^)- 

Proof. First, let us observe that 7\/2 is well defined, as for every / E 36~ ll2 [D.) 
there exists a unique v solution of 12.11) , moreover T\/ 2 is evidently linear. If 
v{x,0) = T\i 2 (f), where v is the solution of (27TTJ, then L\uv = d v v{x,0) and 
from ( |2. 1 1) , L\i2T\i2{f) = d v v{x,0) = fix), i.e. Ti/ 2 is the inverse of the operator 

Il/2- 

In order to show that T\/ 2 is compact when restricted to J£ 2 {Q), let us take f n e 
i? 2 (Q) weakly converging to / £ £6 2 {D.) and consider 7i/2(/n) = v n [x,0) with 
v n e 36 1 { c £) sequence of solutions of JzTTJ with datum /„. From the weak for- 
mulation of l |2.11| l we obtain that v n is uniformly bounded in 36 1 (^) , so that it 
weakly converges to a function v e 36 1 which turns out to be a weak solu- 
tion with datum /. Choosing as test function y/= v n -v\w the equation satisfied 
by v n and taking advantage of the compact embedding of 36 1/2 (fi) in i£ 2 [D.) im- 
mediately gives the strong convergence of v n to v in 36 Y Cio) . And by continuity 
of the trace operator, v n {x,0) = Tmifn) converges to v[x,0) = T\i 2 {f) in ££ 2 (Q). 
Arguing as in Proposition 2.12 in [7] it is easy to obtain that Ti/2 restricted to 
5£ 2 {V£) is self-adjoint and positive. 

Finally, the last part of the statement can be proved by following the argument 



of Proposition 2.12 in |7| (see also Remark 2.13). □ 



To end this section, we face some regularity issues. As already observed, any of 
the above harmonic extensions is of course smooth inside c €. On the other hand, 
improved regularity up to the boundary seems to be prevented by the fact that 
is only Lipschitz. Nonetheless, we can exploit the homogeneous Neumann 
condition (together with some regularity of 3Q) in order to suitably extend the 
harmonic functions outside c €, thus removing that obstruction. 

Proposition 2.19. LetCl be of class C 2 ' a , andlet f £ 36~ ll2 {D.), ve36 l { c £) satisfy 
<2Ju) . Then v £ C 2 ' a [D. x (0, +00)) and 

(i) iff e IP (12), 2 < p < 00, then v £ W^Ptf) and || v\\ w i, P < C(Q, p) II /II lp; 

(ii) iff e W l 'P{n), 2 < p < 00, then v £ W 2 'P{^) and || v\\ W 2, P < C(Q, p)\\f\\ w i, P ; 

(iii) iff £ C°' a (D) then v £ C 1,a (¥) and \\ v\\ c i.a < C(Q,a)||/|| c o,a. 
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Proof. The fact that v e C 2,a (Q x (0, +00)) follows from standard regularity theory 
for the Laplace equation with homogeneous Neumann boundary conditions on 
smooth domains. As far as (i) is concerned, due to the exponential decay of v 



given by Lemma 2.15 we are left to prove regularity near {y = 0}. To start with, 



for any xo e Q, let us consider any half-ball 

5 + = B+(x ,0) = {(x,y): |(x- x ,y)| < R, y > 0} c <€ 
and let us introduce the notation 

H^°iB + ):={^eH\B + ):ifr\ dB+n{y>0] = 0}, a(v,y/) := f Vv-Vy/dxdy. 
Since v solves 12.10) , integration by parts yields 
(2.15) a(v,y/)= I f{x)y/{x,0)dx = I /(x)i//(x,y) da 

JB+n{y=0)} ' JdB+ 

= -J f{x)d y y{x,y)dxdy, Vi// e H 1;0 (5 + ). 

As a consequence, Theorem 3.14 in |27| implies that v e W i,p {B^){xo,0) for e- 
very r < R. On the other hand, let xo e dQ. By assumption, there exists an open 
neighborhood U 3 (xo, 0) and a C 2,a -diffeomorphism $ between U n {y > 0} and 
By (0, 0) which is the identity on the y-coordinate and such that 4>(x , 0) = (0, 0), 
§>{Un c £) = B + n{X]v < 0}. Let v - yoO -1 . Since v is harmonic we have that v sat- 
isfies an equation like |2.15) on B + n {xn < 0}, where now the bilinear form a has 
C l,a coefficients (which depend on D. through the first derivatives of $). Accord- 
ingly, the conormal derivative of v on {xjv = 0} vanishes. Since D. is C 2 ' a , the last 
fact allows to extend v to the whole B + by (conormal) reflection, at least when 
the initial neighborhood U is sufficiently small; in a standard way, the extended 
function satisfies again an equation like (2.15) , and now the corresponding a has 
Lipschitz-continuous coefficients. Furthermore, the analogous extension of / is 
again LP . As a consequence, Theorem 3.14 in |27) implies also in this situation 
that v, and hence v, is W 1,p {B^) for r < R. Taking into account the previous dis- 
cussion, property (i) follows by a covering argument. Finally, (ii) and (iii) can be 
proved with minor changes in the previous argument, by using Theorems 3.15, 
3.12 and 1.17 in (13. □ 

The previous proposition implies a number of regularity properties for the in- 
verse operator 7i/2- Analogous arguments yield improved regularity also for the 
direct operator Z4/2. 

Proposition 2.20. LetD. be of class C 2,a , and let «e jif 1 12 {Q.) besuch that 

ueC l ' a {D.) and d v u = on 3D.. 
Finally, let v e 36 & ) be the Neumann harmonic extension of u according to 



Lemma 2.5 Then v e C 1 '"^), and || v\\ c i,ar^ < C(t2, a) || a|| c i, t 
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Proof. It is sufficient to show that w{x,y) := v{x,y) - u{x) is C l ' a i^€). This can 

once one 



be done by following straightforwardly the proof of Proposition 2.19 
notices that, instead of equation 1 12.15) , w satisfies w[x,0) = and, as v is har- 
monic, 

a(w,iff) = I -S/ x u{x)-'S/ x if/{x,y)dxdy Vi//e Hq{B + ), 
Jb+ 

for every B + c c €. Hence the role that / had in the aforementioned proposition 
is now played by V x u. Since V x u is C '"^), the proposition follows again by 
applying 1 27 1, Theorem 3.12, to w (or to suitable extensions w near 3D. x {0}). □ 

As a conclusion of this section we state the following result, which will be useful 
in the applications 

Corollary 2.2 1 . Let us define the spaces 
X:= {we C 1,a: (Q) : A n u = 0, d v u{x)=0on dn\,F:= j/e C°' a (Q) : A n / = ol. 
Then the operators 

L 1/2 :X^F, T V2 :F^X 
are linear and continuous andL\/2 ° T\n - T\ii°L\ii = Id. 



Proof. The conclusion easily follows from Propositions 2.18 and 2.19 



In the following we will be concerned with positive solutions of equations in- 
volving the fractional operators defined above. In this perspective, the argu- 
ments we employed to improve regularity allow to check the validity of suitable 
maximum principles and Hopf lemma. In particular, the following strong maxi- 
mum principle holds. 



Proposition 2.22. Let c e L°°(Q) and nonnegative. Every u e C l ' a (Q) satisfying 
(2.16) 



'(-A^) 1/2 u + c{x)u> ind 
u>0 in n, 



is either identically zero or strictly positive on O. 

Proof. Let us write u = Z^u + A&u =: u + c u , and let v denote the Neumann 
harmonic extension of u to c €. Then v{x,y) := v{x,y) + c u is harmonic, non- 
negative, and v{x,0) = u{x). Now, if u[xq) = for some Xq e Q, this would imply 
d v v{xo,0) = (-A) 1/2 m(xo) > 0, in contradiction with the Hopf principle for har- 
monic functions. On the other hand, if u{xq) = for some Xo £ d£l, we can argue 
in the same way, considering instead of v its conormal even extension, as in the 
proof of Proposition|2.19| □ 



Remark 2.23. As a direct consequence of Proposition 2.22 and of classical max- 
imum principle for harmonic functions, we deduce that, if u > satisfies l |2.16) , 
then its harmonic extension v- v + c u is positive in c €. 
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3 The Weighted Logistic Equation 



Our main application is the study of the positive solutions of the nonlinear Prob- 
lem {L7j, understood in terms of the operator (-A j/) 112 . To this aim, a necessary 
solvability condition is that the right hand side of the equation has null average. 
On the other hand the possible solution u, being positive, has positive average. 
In order to apply the theory developed in the previous section, we recall that any 
u e i? 1/2 (Q) can be decomposed as 



(3.1) 



u = An u + Zf l u=: u + c u 



where c u is constant and u e 36 ' (CI). Using Lemma 2.5 we can denote by v the 
Neumann harmonic extension of u to , obtaining that 

v{x,y) := v{x,y) + c u 

is harmonic and v{x,0) = u(x). It is worthwhile noticing that, as far as c u + 0, 
vt H 1 ^). 

Taking into account the previous discussion, we can now define what we mean 
by a weak solution of a general nonlinear problem. 



Definition 3.1. A weak solution u: Q 
(3.2) 



: of the nonlinear problem 



(-A) 1/2 k = /(x,u) inO, 
d v u = on 3D, 



is a function u e H ll2 {Cl) such that /(•, u(-)) e H" 1/Z (Q) and 



112, 



both A n f[-,u) = and (-A^) 1/2 k = /(•, u). 

In particular u{x) = v (x,0), where v{x,y) = v{x,y) + h, and {v,h)£ J€ l ^€) x U is 
a weak solution of the nonlinear problem 



(3.3) 



in the sense that 



Av = 
d v v = 

d v v- f{x, v + h) 



in^, 

on dQ x (0, +oo), 
on Q x {0}, 



f 

Jo 



fix, vix, 0) + h)dx = 0, 



(3.4) 



/ VvVif/dxdy= \ fix,v+h)y/dx Vif/E^i'ig), 

I fix, vix, ( 
Jn 



,0) + h)dx = 0. 
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Using the previous definition we can now rewrite Problem (L7) in the equivalent 
form 



(3.5) 



(-A jf) ll2 u = Xu(m(x) - u), 



X I u{m{x) -u)dx = 0, 
Jn 

and we recall that we assume 

X>0 and meC ' 1 ^). 

Remark 3.2. In the standard diffusion case, nonlinear boundary data have been 
frequently considered especially in the determination of selection-migration prob- 
lem for alleles in a region, admitting flow of genes throughout the boundary (see 
fTBl and the references therein). 

As in the classical literature concerning the logistic equation, the comprehen- 
sion of the linearized problem arises as crucial in the study. In our context, this 
correspond to tackle the following weighted eigenvalue problem 



(3.6) 



(-A j/) ll2 u = Xm{x)u, 

X I m(x)udx = 0. 
Jn 



Remark 3.3. When m=l, the nontrivial solutions of 

{-A^) V2 cp = Xcp, 

I (pdx = 0, 
Jn 

are (Pk = 4>k associated to = y/pfc where (p^ and jj.^ are respectively eigen- 
functions and eigenvalues of the usual Laplace operator - A with homogeneous 
Neumann boundary conditions as in (1.2) . 



Remark 3.4. Taking into account the usual decomposition u = u + c u as in {37T 
we have that Problem 13.6|l can be written as 

u- XTi /2 (m{u+ c u )) = 0, 

/ m{u + c u ) dx = 0. 
Jn 



where T\/2 is compact by Proposition 2.18 If moreover we assume f n m^0, we 



can solve the second equation for c u and infer the equivalent formulation 



u — XTi 
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- fn ma 
mu m 

fn m 



0. 



As a consequence, we can apply Fredholm's Alternative, obtaining that the spec- 
trum of the operator at the left hand side consists in a sequence of eigenvalues 
(Ajfc)fc, with associated kernel of dimension d^ < +oo and closed range having 
codimension d^. 
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Lemma 3.5. Any nontrivial solution u of Problem is of class C l ' a {D.) . More- 
over, u>0 implies u> in Q. 

Proof. The proof relies on the classical bootstrap technique. Indeed, as above, 
let us write u = u+ c u and let us denote with v the Neumann harmonic extension 
of u. Since m e C°' l {D.), Proposition 



2.19 



and the trace and Sobolev embedding 



theorems imply 



ue W 1 " 1/r ' r (0) 
==> Amu e L 



ue L 

Nrl(N+\-r) 



Nrl(N+l-r) 



(Q) =^> DeW 



l,Nrl(.N+l-r) 



whenever 2 < r < N+ 1. Starting from ro = 2 and iterating the above procedure 
the first part of the proposition follows. As a consequence, the second one is 
implied by Proposition|2.22| □ 



Searching for positive solutions of (3.5) , we are interested in positive eigenfunc- 
tions of \3.6\ . Of course, A = is always eigenvalue with normalized eigenfunc- 
tion (po = H \ft&l > 0, but this does not prevent the existence of positive eigen- 
functions associated with positive eigenvalues. 

Lemma 3.6. If there exists a positive eigenvalue Ai with a positive eigenfunction 
(pi then the function m is such that 



(3.7) 



/ m(x)dx < and 3 xo e n such thatm{xo) > 0. 
Jn 



Proof. Supposing that there exists Ai > with positive nonconstant eigenfunc- 
tion (p\,we can apply Lemma : 



and use cpo = l/ vl^I > as a test function in 



2.15 



the weak formulation of \3.6\ satisfied by (p\ to obtain 

= A 1 | m{x)q)i{x) dx. 
Jn 

As <f>\ is positive, m has to change sign. Now, taking advantage of the usual de- 



composition, let us write cpiix) = V\ {x,0) = Vi{x,0) + h. From Remark [2 .23 
deduce that V\ > on S^, and Lemma 



2.15 



we 



allows to use y = 1 / v i as test function 



in the equation satisfied by <f>\. We obtain 

2 



i 



Mvi 



dxdy = X\ I m{x)dx, 
Jn 



and the lemma follows. 



The following result shows that the previous necessary condition is also suffi- 
cient in order to obtain the existence of a first positive eigenvalue with positive 
eigenfunction. 



Theorem 3.7. Let us suppose that m e C • (Q) satisfies condition (377). Then 
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(i) there exists Ai > and <p\ e C 1,a (Q), <p\ > in Q, solution of with 
A = X\; 

(ii) Ai is the smallest positive eigenvalue, it is simple and any other solution of 
(3 .6) with A > Ai changes sign. 

Proof. We will find cpi solving the extension Problem l |3.3| l, via a constrained 
minimization. Namely, we look for a minimum of the functional E: 3€ x &>) x 
U — U defined by 

E{v,h) = - f \Vv\ 2 dxdy 
2 JV 

constrained on the manifold 

(3.8) M = {{v,h)eje l { c £)xU: [ m{x){V{x,0) + hfdx = l}. 

Jn 

First, let us observe that M + 0, indeed, from 0.7\ we can find a) an open set 
of positive measure such that m{x) > in to, so that for any w e C^°(oi) having 
zero mean there exists a suitable positive real number C, such that {wl \/C,0) 
belongs to Jt. Let us consider a sequence {v n , h n ) e Jl such that 

E[v n ,h n ) ^inf£>0. 

M 

From the definition of E it follows that v n is uniformly bounded in J6' l { c £), so 
that v n {x,0) is uniformly bounded in H ll2 {£l); by the compact embedding of the 
trace space H ll2 {£l) in L 2 (D) we obtain that there exists V\ e JC 1 ^) such that 
v n tends to v\ weakly in J? 1 ^), and v n {x,Q) strongly converges to z>i(x,0) in 
L 2 {D.). As far as the sequence h n is concerned, let us show that it is bounded by 
contradiction, assuming that, up to a subsequence, h n -* +oo (the case h n -* 
-oo can be handled analogously). By the definition of M it follows 

lim h n \ m{x)dx = -2 \ m{x) v n (x,0)dx + o(l), 

n^+oo J n Jn ' 

where o(l) denotes a quantity tending to zero as n goes to infinity. Then, a con- 
tradiction follows from l |3.7) . As a consequence, there exists hi such that h n -> h\ 
in R. By weak lower semicontinuity of E, it results that the pair {v\, h\) satisfies 

(3.9) E{vi,h\) = infE, / m(x)(i5i(x,0) + hifdx = 1. 

^ Jn 

Moreover, let us show E{V\,hi) > 0. Again by contradiction, let us assume that 
E{v\, hi) = 0; then, as V\ £ ^^f 1 ^), it follows £>i = 0, and from l |3.9| l we obtain 
that 

h\ \ m{x)dx = 1, 
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so that \3.7) yields again a contradiction. Since {v\,hi) is a constrained mini- 



mum point of E on Ji, there exists A e R such that, by Proposition 2.15 for every 
y/ e ^(^n satisfying (2A2), it holds 

/ V£»i(x, y)Vi//(x, y)dxdy = A / m(x)(£>i(x,0) + fti)i//-(x)<ix, 
J^f JO 



(3.10) 
(3.11) 

Choosing y/ = v\ + h\ we infer 



/ m(x)(i>i(x,0) + h\)dx = 0. 



inf£ = E{v\,hi) = \ \Vi>i\ 2 dxdy = A / m(x)(z)i(x,0) + fai) 2 <ix = 



A, 



thus A > and we can define 

(3.12) Ai := A = inf£ = E{v, hi) and </>i(x) := v\{x,G) = v\[x,0) + hi 



its corresponding eigen func tion, which is a weak solution of Problem l |3.6[ . As 
a consequence, Lemma 3.5 implies that <p\ e C 1,a (Q) for every a e (0, 1). Since 



implies that q>\ 

any other solution (A, v) with A > of (3.6| l corresponds to a constrained crit- 
ical point of E on M, X\ is the smallest positive eigenvalue. In order to show 
that ip\ can be chosen positive, let us take w{x) = \(p\[x)\ = |£>i(x,0) + h\. Writing 
w{x) = w{x) + c w , with w e J¥? ll2 [Q.) and c w constant, let us consider £(x, y) e 
the harmonic extension of w{x) obtained thanks to Lemma 2J5 No- 
tice that {(,c w ) e M; moreover ({x,0) = \ i>i[x,0) + h\ — c w = z(x,0) for z(x,y) = 



I V\ (x, y) + h\ - c w so that, thanks to Remark 2.7 



E(w)<E{z) 



= ( |V| 



vi{x, y) + h\\ 



■■ ( IV 



v l \ 2 = E{v l ). 



As a consequence, also the nonnegative function w solves the minimization 
Problem (3.9| l, showing that we c an a ssume, without loss of generality, that q>\ 
is nonnegative. But then Lemma 



3.5 



applies again, yielding cpi > on Q. It is 
possible to show that Ai is simple by contradiction, supposing that there exists 
(pi and u solutions of J3.12) , with </?i(x) = yi(x,0) = £i(x,0) + h, u{x) = w{x,0) = 



w{x, 0) + k. From Remark 2.23 we deduce that v\ (x, y) > in , so that we can 
use if/[x,y) = w 2 (x, y) I V\ (x, y) as test function in the equation satisfied by cp\, 
obtaining 



Ai / m(x)u 
Jn 

This implies 



[x)dx= f V^i(x,y) 



w{x, v) 
2 - \ Vw{x,y)- 
vi (x, y) 



' w[x,y) 
v\{x,y)l 



Vviix, y) 



dxdy. 



Ai / m{x)vf{_x)dx - 
Jn 



f I w{x,y) 2 
-/ \\7w{x,y) Vyi(x,y) dxdy 



fi(x,y) 



+ Ai / m{x)u' L {x]dx, 
Jn 
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that is 



= - / |v 



w{x, y) ■ 



w{x,y) 
vi {x, y) 



Vfi(x,y) 



dxdy 



J<€ 



y) 



. ' w{x,y) ' 



dxdy, 



yielding the linear dependence between cp\ and u. Moreover, it is possible to fol- 
low the same argument as in 1 19 1 to obtain that also the algebraic multiplicity of 
Ai is one. 

Now we come to part (ii). In order to show that there is not a positive solution 
u of Problem (3 .6) associated to A > X\, let us argue again by contradiction, and 
suppose that there exists u[x) = w (x,0) = w[x,0) + k positive eigenfunction as- 



sociated to an eigenvalue A greater than Aj. As before, observe that Remark 2.23 
allows to choose as test function y/{x, y) = V\ (x, y) 2 lw{x, y) in the equation sat- 
isfied by u and obtain 



A / m{x)(p\{x)dx= \ Vw{x,y) 
Jo. 

which gives 

(A-Ai) I m{x)ip\{x)dx 
Jn 



vi (x, y) 
2 1 \ Vviix,y)- 
w{x,y) 



' v\{x,y) 
w{x, y) 



Vw{x, y) 



dxdy, 



= -/|v 



vi (x, y) • 



vi(x, y) 
w{x, y) 



Vw{x, y) 



dxdy 



and j3.9| l immediately implies that A < Ai. 



□ 



Remark 3.8. Notice that changing m in -m we can prove that there exists a first 
eigenvalue A_i < with a positive eigenfunction <p\ when the following condi- 
tion holds 



(3.13) 



/ 

Jn 



m{x)dx>0, 3x e Q, such that m(xo)<0. 



Moreover, arguing as in the end of Theorem 3.7 it is possible to prove that there 
are not positive eigenvalues with positive eigenfunctions, when m has positive 
mean. 



We are finally in the position to tackle the logistic equation l |3.5) ; let us start our 
study with the easy observation concerning the autonomous problem, i.e. m = 
1, contained in the following proposition. 



Proposition 3.9. If m = 1 then every non-negative solution u of \3.5) is either 
u = or u = 1. 

Proof. Let u[x) = v{x,0) = £*(x,0) + c u be a solution of J33). Lemma [2.15| implies 
that we can choose as test function y/{x, y) = v[x, y)- 1. Since Vifr = Vvwe obtain 

0< / \Vv{x,y)\ 2 dxdy = -A / u{x){u{x) - ifdx < 0, 
J<€ Jn 



and the lemma follows. 
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Back to the nonautonomous case, the next two results provide a priori bounds 
on the set of the positive solutions of (33) and on the set of the parameters X. 

Lemma 3.10. Let 

M = sup m + . 

Q 



Then any positive solution of (33) satisfies 

(3.14) u{x)<M. 

In particular, ifm[x) < then no positive solution exists. 

Proof. Let u{x) = v[x,0) = v{x,0) + c u be a positive solution of (33). We use 
Lemma 



2.15 



to choose as test function q){x, y) = {v{x, y) - M) + and obtain 



\ \Vv{x,y)\ 2 = X \ f{x, u(x)){u(x) - M) + < -X [ [(u{x) - M) + fu{x), 

where c & + = {{x,y) e ^ : v[x,y) > M] and/(x, s) = s(m(x)-s). Since the left hand 
side is non-negative, we obtain that [{u[x) - M) + ] 2 u[x) = a.e., and the lemma 
follows. □ 

Corollary 3.11. Any nonnegative weak solution of (33) is C 1,a (0) and strictly 
positive on Q. 



Proof. This is an easy consequence of Lemma 3.10 and Propositions 2.19 and 
IZ221 



□ 



Concerning the set of the parameters X the following necessary condition holds. 



Lemma 3.12. Assume (377). Then, if there exists a positive solution of (33), then 
X > X\. 



Proof. Let u{x) = v[x,0) = v{x,0) + c u , be a solution of equation (33) and let 
V\ = v\ + h satisfying l |3.10| l and (3TTTJ. By Corollary 3.11 we can take as a test 
function in the equation satisfied by V\ , y/ (x, y) = v 2 [x, y) I v\ (x, y) to obtain 



C \ v{x,y) i v{x,y) \2 

VvAx,y) 2—^-Vv 1 (x,y)-[—^-) V Vl [x,y) 
[ vi(x,y) \vi{x,y)> 



v\ [x, y) 
and by using (33 



dxdy = Xi \ m{x)u 2 {x)dx 



> — Ai / u 3 {x)dx- I \Vv(x, 



v{x,y) 
v\{x,y) 



Vvxix, y) 



dxdy 



X\ 



[ y\v{x, 



y)\ dxdy, 



providing the conclusion. 
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We will obtain existence results for Problem {33) via Bifurcation Theory; deve- 
loping this approach we have to take into account that every solution may have 
a constant component that is invisible in the differential part of the equation, 
then in order to make this component appear, we will be concerned with the 
map G:IRx[RxX— ► Y where X is defined in Corollary |2.21| Y is defined as 

Y = e C°' a (Q) xU, t = J w{x)dx,\ 

and G has components Gi(A, ft, u) and G2(A, ft, u)) given by 

Gi(A, ft, u) = -L\i2{u) + Af[x,u + ft), 

G2(A, ft, u) = f G\{X,h,u)dx= \ Af{x,u+h)dx, 
Jn ' Jn 



(3.15) 



for f{x,s) = s{m{x) - s). Let us remark that, since f a G\dx = G2, we have that the 
elements in the range of G automatically satisfy the condition in the definition 
of Y. Moreover, thanks to Corollary|3.11| the zeroes of G correspond to solutions 



of Problem |3]5). Of course, we are interested in nontrivial solutions. 
Definition 3.13. We denote the sets of trivial solutions of G(A, ft, u) = as 

&\ := {(A, 0, 0) : A > 0} , ST 2 := {(0, h, 0) : h > 0} , 
and the set of positive solutions as 

S?:= Ux,h, u) : A>0, w + ft >0 in n} . 

implies SP = 0. 



3.10 



Remark 3.14. We observe that if m < on O then Lemma 
On the other hand, if m is a positive constant, reasoning as in Lemma|3l) 



we in- 



fer that Sf - {(A, m, 0) : A > 0}. As a consequence, in the following we can assume 
without loss of generality that m is not constant and is positive somewhere. 

The following local bifurcation result is concerned with the case of negative 
mean of the function m. 

Proposition 3.15. Let condition \3.7\ hold and let Ai be defined as in Theorem 



3. 7 Then (Ai , 0, 0) is a bifurcation point of positive solutions of Problem \3.5\ from 



ET\, and it is the only one. Moreover, locally near such point, £P is a unique C 1 
cartesian curve, parameterized by A e (Ai, X\ + 5),forsome5 > 0. 

Proof. The proof relies on classical results about the local bifurcation from a 
simple eigenvalue, see for example |2 |, Chapter 5, Theorem 4.1. 
The derivative of G with respect to the pair (ft, u) has components 

(3.16) 5(^jj)Gi(A, ft, u)[k, iv] = -L1/2W + A/ S '(x, u + h){w + k), 

(3.17) d(fe,c)G2(A, ft, u)[k, iu] = A / f' s {x,u + h){w + k)dx, 

Jn 
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which, evaluated at the triplet (A, 0,0), gives 



d{h a) G(A, 0, 0) [k, w] = [-L1/2W + Am(x)(u> + k),X I m(x){w + k)dx). 

Jo 



Now, by Remark 3.4 we have that (A, 0, 0) can be a bifurcation point for positive 
solutions only if there exists a pair [k, w) with w + k > belonging to the kernel 
of the operator dy^GfA, 0,0), i.e. such that 

L112W- Xm{x){w + k), 

X\ m{x)(w + k)dx = 0, 
Jo 

which is equivalent to say that the function w{x,y) = w{x,y) + k is a positive 
solution of 



(3.18) 



(-A jy) 1I2 w = Xm{x)w, 

A I m(x)wdx = 0. 
Jo 



For this linear eigenvalue problem, Theorem 3.7 shows that there exists only one 
positive simple eigenvalue Ai with a positive eigenfunction ip\ satisfying l |3.12} . 
Decomposing ip\ as (pi + c Vl = (p\ + hi, we deduce that the kernel of the opera- 



tor d(/j,£)G(Ai, 0,0) is generated by (^1,^1). By virtue of Remark 3.4 this implies 
that the range of the operator d[h,U) G(Ai, 0, 0) is closed and that it has codimen- 
sion one. Such range consists in the pairs {w, t) such that there exists a solution 
z{x, y) = z[x, y) + h of the problem 

(-A ,jy) II2 z = X\m{x)z+ w, 

Ai I m{x)zdx = t. 
Jo 

Taking as test function y/{x, y) = v\{x,y) in the weak formulation of the first 
equation we derive that the range is given by 

I ( w, t) £ Y such that J w {x)(p\{x) dx = o| . 

Deriving \3.16\ and (3.17) with respect to A leads to 

dxd(h a)G{X, h, u)[l, k, iu] = [f' s {x,u + h){iv + k), \ f' s {x,u + h)(w + k)dx) 

Jo 

and, denoting with M the operator 6x3^,0} G{X, h, u) [A 1( 0, 0] we have that 

M(hi,<pi) = [m{x)[(p\ + hi), \ m{x)((f>i + hi)dx) = [m{x)(pi, I m(x)cp\dx). 

Jo Jo 
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At this point, in order to apply the aforementioned theorem from |2|, we only 
have to check that M{h\,ipi) does not belong to the range of d^,a)G, and this 
occurs because 

( m{x){(pi{x) + h\){.ipi{x) + h{]dx = [ m{x)ip\dx = — [ |V</>i| 2 = l. 
Jn Jn M J<e 

Then at (Ai,0,0) a bifurcation occurs. Moreover, as f{x,s) = s[m[x) - s) is of 
class C 2 with respect to s, the set of the nontrivial solution of G(A, u, h) = near 
(Ai,0,0) is a unique C 1 cartesian curve, parameterized by 

(3.19) A = Ai + //(?), h = thi + j6(Ai + p,{f), t(p\), u= t(pi+y{Xi + fj,{t),t(pi) 

for t e (-£,£), t # 0. Here both + /i(J), ?^i) and jS(Ai + ^x(f), f^i) are o[t) 
as ? — ► 0, while a direct computation shows that ji/(0) > 0. Thus, for sufficiently 
small t > 0, it is possible to write t = t (A), and the solution (A, h, u) is positive. □ 

Coming to the case of positive mean of the function m, it is more convenient to 
use as a bifurcation parameter h instead of A. 

Proposition 3.16. Assume 

(3.20) | m[x)dx>0, 



and let h* be defined as 



Jp. 

h* = l 



(3.21) h* = f m{x)dx. 

Jn 

Then (0, h* ,0) is a bifurcation point of positive solutions of Problem \3.5\ from 
&2, and it is the only one. Moreover, locally near such point, Sf is a unique C 1 
cartesian curve, parameterized by A e (0, S] , for some 5 > 0. 

Proof. The derivative of G with respect to (A, u) has components 

(3.22) 5(^ a) Gi(A, h, u)[l, w] = -L1/2W + A/j(x, u + h)iv+ lf{x,u + h) 

(3.23) d(x s)G2(A, h, u)[l, w] = A / f' s {x,u + h)wdx+l I f{x,u + h)dx, 

Jn Jn 

so that a pair [l, w) belongs to the kernel of fly,,^) G(0, h,0) if and only if (/, w) 
solves the problem 



(3.24) 



-L1/2 w + Ifix, h) = 0, 

I [ f{x,h)dx=0. 
Jn 



For / = 0, taking into account Corollary 2.21 we find w = 0, while for I + the 
mean of f{x, h) has to be zero and this, thanks to (3.20| l, yields the positive value 
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for h* given by (3.21) . With this choice of h* , for any I there exists a unique so- 
lution of the first equation. Denoting with z* the one corresponding to I = 1, we 
obtain that the kernel of d^,a)G{0, h,0) is the one dimensional space generated 
by the pair {1,2*). 

On the other hand, pair ( w, t) belongs to the range of d(x,U)G(Q, h* , 0) if and only 
if there exists a solution [v, I) of the problem 

Li/2V= lf{x,h*) + w, 

I [ f(x,h*)dx=t. 
Jn 

Since the function f[x, h* ) has zero mean, t has to be zero and the range is given 

by the set {{w, t) e Y, such that t = 0} which is closed and of codimension one. 

Deriving (3.22| l and (3.23) with respect to h leads to 

» . 
dfcd(A,fi)Gi(A, h, u)[l, w] = Xf s (x,u+ h)w + lf s (x,u + h), 

dhda «)G 2 (A, fa, u)[l, iv] = A I f s [x,u + h)wdx+l I f s {x,u + h)dx. 

Jn Jn 

This time we obtain the operator N = 0^fl ( A,S)G(O, h*,0), which computed on 
(l,z*) gives 

N{l,z*) = {f!.{x,h*), f fl{x,h*)dx) 

If the second component of N{l,z*) + then N(l,z*) does not belong to the 
range of ,n)G{0, h* ,0) implying that bifurcation occurs in this case too; and 
this is true since (3.21) yields 

N{l,z*) = \m(x)-2^ m{x)dx, - j m{x)dx^. 

As before, using |2 |, Chapter 5, Theorem 4.1, we deduce the existence of a carte- 
sian curve in a neighborhood of (0, h* , 0) with representation 

h=h* + v{t), X=t + a{h* +v{t),tz*), u= tz* +j{h* +v{t),tz*), 

for t £ {_-£,£), t±0. Here both a{h* + v{t], tz*) and j{h* + v{t), tz*) are off) as 
f -» 0, while v(0) = 0. Since h* is positive and also A is positive for t positive and 
small, the proposition easily follows. □ 



Remark3.17. We stress the fact that both in Proposition 3.15 and in Proposition 



3.16 we can locally parameterize Sf with respect to A, even though in the latter 



the bifurcation parameter is h. 



Remark 3.18. If (3.13) holds we can go through the proof of Proposition 3.15 and 
use Remark 3.8 to obtain that X\ < is bifurcation point of positive solutions of 
(3.5| l with A < and u = u+ h nonnegative. Moreover, as in the case of Ai > 0, 
Lemma 3. 12| implies that the bifurcation occurs on the right hand side of Ai. 
Finally, let us notice that, in order to show the local bifurcation from (0, h* , 0), it 
is enough to assume (3.20) and m needs not to be sign-changing. 
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Note that, by Proposition 2.18 it is possible to reformulate the equation G = 
in terms of a identity minus compact map, see also Remark 3.4 Then a classi- 
cal result due to Rabinowitz |22| implies that the continuum bifurcating either 
from (Ai,0,0) or from (0, h* ,0) is actually global. Here we prefer to recover this 
result from a stronger one: indeed we are going to show that the set £f of positive 
solutions is a smooth arc. 

Lemma 3.19. Let (Ao, Uq, ho) e SP. Then there exists e R x U x X neighborhood 
o/(Ao, uo, ho), 5 > and a C 1 map W : (Ao - 6, Ao + 5) —- U x X such that 

Sf> n °U = {(A, (A), W 2 (A)) : A e (A - <5, A + 5)} . 

Proof. The conclusion will follow from the application of the Implicit Function 
Theorem to the map G(A, h, u) defined in J3A5). To this aim, taking into account 
13.16| l, J3A7), we want to show the invertibility of the operator 



-L1/2Z- Ao(m-2uo)(z+ t),Ao J {m-2uo){z+ f)j . 



d(h,u)G(h,ho, u )[t, z] 

We claim that, for such operator, the Fredholm Alternative holds. Reasoning as 
in Remark 3.4 to obtain the claim it is enough to show that f n {m{x) - 2uq) t 0. 
But this can be easily obtained by testing the equation for uo with 1 / vo, where as 
usual uq (x) = uo [x) + ho = v (x, 0) + h = v (x, 0) : 



/ 



(m(x) -2uo)dx < \ (m(x) - uo)dx 



/< 

Jn 



-f\- 



vo 
vo 



dxdy. 



Once the Fredholm Alternative is established, we have that d^,u) G(Ao, ho, Uo) is 
invertible if and only if its kernel is trivial. In turn, {t, z) belongs to the kernel if 
and only if z = z + t solves the problem 



(-A^) 1/2 z = A (m(x) -2uo)z, 

Ao I (m(x) - 2uq)z = 0. 
Ja 

Taking y/{x, y) = w 2 (x, y) / Vq (x, y) as test function in the equation satisfied by Vo, 
where w is the harmonic extension of z, we obtain 



Vvoix, y) 



w 

2—Vw{x,y)- 
v 



IV 
\Vo 



Vy (x,y) 



dxdy = Ao I 



(m(x) - Uq[x)) z {x)dx. 



Then we can test the equation for w with w itself, and subtract it from the equa- 
tion above. We obtain 

f I w{x,y) 2 r 

0<-/ \Vw{x,y) Vvo{x,y) dxdy = Aoj uo{x)z {x)dx<0, 

Jn ] ' v (x,y) ' Jn 



which implies that z, and then (?, z), must vanish. 
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Theorem 3.20. Let Sf he as in Definition 3.13 Then 



(i) if H3.7| l holds then &> is the graph of a C 1 map Y : {X\, +00) — R x X, with 
¥(A+) = (0,0); 

(ii) if ( |3.20) hoZds then SP is the graph of a C 1 map *P : (0, +00) -^IxX, with 
Y(0 + ) = (fc*,0). 

Proof. To start with, we prove that ^ contains such a graph. Let us assume con- 
dition J3.7| l, and let us define 

A:= sup{3^£C 1 ((A 1 ,A),RxX), graphOP) c Sf\ . 

A>A, 



Proposition 3.15 and Lemma 3.12 imply that A > X\, let us suppose by contra- 
diction that A < +00, and consider a cartesian curve W: (Ai,A) -»RxI, de- 
fined by ¥(A) = ( v Pi(A),¥ 2 (A)) ) with {X,W 1 (X),V 2 (X)) e 5?. Let us consider a se- 
quence X n < A tending to A with corresponding solutions {X n ,h„, u n ), where 
h n = ^1 (X n ), u n - X I / 2(A«), and u n = u n + h n . Moreover, let us recall that we can 
write u n (x) = v n {x,0) and v n (x,y) = v n {x,y) + h n . Taking as test function in (3.5) 
ipix, y) 
bound 



v n [x, y) and applying Lemma 3.10 we immediately infer the uniform 



(3.25) 



from which we deduce that u n is uniformly bounded in the spaces L p [D.) with 
1 < p < 2N/{N— 1). Since G2(X n , h n , u n ) = 0, where G 2 is defined in l |3.15) , and as 
u n has zero mean, h n has to be positive and we obtain 



h 2 n \n 



< h n / m{x)dx-2\ u n (_x)dx 
Un Jo. 



+ L<cL+ h n M\Q\, 



for c positive constant and M defined in Lemma 3. 10 Hence, also h n is bounded 
and there exists h > such that, up to subsequences, h n -» h, and u n = u n + 
h n — - u + h > 0. From Proposition 2.22 we have two possibilities, either u > or 
u = 0. In the first case we have obtained a positive solution of (33) with A = A 



and Lemma 3.19 provides a contradiction with the definition of A. In the second 
case, A turns out to be a local bifurcation point for positive solutions, but then 
Proposition |3 . 1 5| implies that A = Ai which is again a contradiction, showing that 
A = +00. 

When l |3.20) is assumed we define 

A := sup {3 ¥ e C 1 ((0, A), U x X), graph(Y) c S?} . 

A>0 



Then A > by Proposition 3.16 and arguing as above we obtain that also in this 
case A = +00. 

Finally, we are left to show that Sf \ graphs) is empty. We prove it assuming (377) , 
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when 13.20) holds the same conclusion can be obtained with minor changes. 
Let us argue by contradiction and suppose that there exists A* with distinct 
positive solutions {A* ,h\,u\) and (A*,/j2,U2)- Arguing as above, it is possible 
to see that {A*,hi, U\) and (A* ,h2,U2) belong respectively to global branches 
5P\ and Sfy, of positive solutions that can be parameterized by cartesian curves 
Wi, W 2 : [Ai, +00) — U x X. Notice that ^ n ^2 = and neither SP\ nor may 



have turning points, otherwise Lemma 3.19 would be contradicted. As a con 



sequence Ai is a multiple bifurcation point of positive solutions, but this is in 
contradiction with the local representation provided in < |3. 19) . □ 

Corollary 3.21. There exists exactly one positive solution u{x) = u{x) + h associa- 
ted to any A > Ai when \?>.7) holds, and there exists exactly one positive solution 
u + h for every A > when \3.20\ holds. 



Proof. This is an evident consequence of Theorem 3.20 



Taking into account Remark 3.14 we have that the only case left uncover by The- 
orem 3.20 is when m has zero mean but it is not identically zero. Notice that 
in such case the candidate bifurcation point is the origin, but it is not possi- 
ble to argue as in the previous results, as the mixed derivatives dxd(h,tt)G(Q,0,0), 
^/i^(A,s)G(0,0,0) are now both trivial. Nevertheless, we can still prove the exis- 
tence of a solution for every A > arguing by approximation. 



Theorem 3.22. Assume that m is a Lipschitz function not identically zero and 
satisfying 

\ m{x)dx = 0. 
Jo. 

Then SP is the graph of a C 1 map W : (0, +00) -»RxI, with W{0 + ) = (0, 0). 



Proof. Let us choose no > 1 such that for every n > no the weight 

1 

m n [x) = m(x) 

n 

satisfies hypothesis 13.7| l. Letting Jl n be defined as in 13.8| l, with m n in the place 
of m, Theorem |3.7| yields the existence of a first positive eigenvalue 

X\ n - inf / \Vv(x,y)\ 2 dxdy 

associated to the weight m n {x). Let us define 

9 

u n {x):=p n m{x) + q n , where p n := 



Vnf n m 2 ' 

We claim that, when n is sufficiently large, q n can be chosen in such a way that 
u n s.M n , that is, 



J. 

Jn 



m„{x)ip n m{x) + q„fdx = 1. 
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Indeed, since f n m 2 = f a mm n , by direct calculations the above equation can be 
rewritten as 

, ^_ 4f C) m 2 m n 
q 2 n -4s/iiq n + n- ff = 0, 

which solvability is equivalent to the condition 

[fnm 2 f 

trivially satisfied for n large. With this choice of q n we have that, denoting with 
rh the Neumann harmonic extension of m and writing v n [x, y) = p n m{x, y) + q n < 
it holds v n {x,0) = u n {x). This implies 

X\ n < I \Vv„(x,y)\ 2 dxdy = p 2 n \ \Vm(x,y)\ 2 dxdy, 

yielding X\ >n — as n —- oo. 
Now, Theorem 



3.20 



i provides a sequence of C 1 functions : [Ai, n , +oo) — ► R x X 
with ^(A) = [h n ,u n ) positive solution of (3. 5) with weight m n {x). Let us fix 
< 5 < A and n\ > no such that X\ iH < 8 for every n > n\, so that W n is defined in 



[8, A] for every n > n\. Using Lemma 3.10 and Proposition |2.19| we obtain that 



{h n , u n ) is uniformly bounded in K x X, so that up to a subsequence {h n , u n ) con- 
verges in U x 1 Cg) to a pair {h, u) solution of (3.5) ; moreover, the same a priori 
bounds implies that \P „ satisfies the hypotheses of Ascoli-Arzela Theorem in the 
closed, bounded interval [8, A] , yielding the existence of a continuous function 
*P : [8, A] such that ¥ n converges to W uniformly and ^(A) = {h, it). By the arbi- 
trariness of 8 and A, we have that W is defined in the whole interval [0, +oo), and 
T(0) = (0,0). 

The only thing left to show is that ^(A) + for ever A > 0. Let us argue by contra- 
diction and suppose that there exists A > such that ^(A) = {h n , ii n ) -> (0,0). 
As usual, let u n = u n + h n and v n = v n + h n be such that v n {x,0) = u n (x). Setting 
z n = u n /\\ v n \\^i m , w n = v n l\\ VnWjpucg), we obtain 

/ Vw„Vy/ = xl z n (m„ix)-u n )y/, I z„(m„(x) - u n ) = 0, 
J<r? Jn Jn 

for every test function y/. Passing to the limit we obtain 

/ VwVy/ = Xf m(x)zif/, I m(x)z = 0, 
J^ Jn Jn 

which is equivalent to say that the nontrivial function z is a nonnegative eigen- 
function associated to the positive eigenvalue A, but as m has zero mean value 
this contradicts Lemmas|3.5||3.6[ □ 
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Figure 1: graphs of A(s, m\, (0,L)) (dots) and A(s, m2, (0,L)) (squares) as functions 
of s £ [0.4, 1], for L = 2.5 (top left), 1 = 5 (top right), 1 = 8 (bottom left). At bottom 
right, the graph of A{s, m\, (0,1)) for L = 3.5. 



As we mentioned in the introduction, a relevant question is the one of compar- 
ing the two eigenvalues 

Ai(m,Q) = inf-j I \Vv{x,y)\ 2 dxdy: I m{x) v 2 {x,0)dx = 1, / m{x)v{x,0)dx = 0> , 

Jiy(m,D.) = inf < \ \Vu{x)\ 2 dx: \ m(x)u 2 {x)dx = 1, / m{x)u[_x)dx - >, 
[Jn Jn Jn J 



which correspond to the linearized version of (L7J and \l.8\ , respectively. We 
conclude this section showing some simple numerical results in dimension N = 
1. Using the usual Fourier representation with basis defined in \1.2) , we have 
that 



mu 



=1 



J m{x) |^ Mj(^y(x)| <pi{x)dx 

X M j / m{x)(pj{x)(pi[x)dx 
j \Jo , 



L M ij u i 

j 



Under this point of view, solving the above minimization problems amounts to 
finding the smallest positive eigenvalue A(s, m, O) of the problem 

diag(MD/>o u = AMu - 

indeed \\{m,£l) = A(l/2, m,Q) and ^(m,!!) = A(l, m,Q). In turn, such eigen- 
value can be easily approximated by truncating the Fourier series. In Figure[T]we 
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report these approximations in the cases Q = (0, 1) and 



m\(_x) = cos 




(271 ) 1 

m?{x) = cos — x — . 

\L } 2 



Hence m has always mean equal to 1/2, while 



71 

Ml (L) = -2. 



We observe that in the case L = 2.5 < n then /^i > 1, and thus A is increasing in 5 
for any choice of m as one can trivially prove. On the other hand, when /ii < 1 the 
situation is more variegated. In any case, the eigenvalue corresponding to m\ is 
always lower than the one corresponding to wi2, in agreement with the results 
obtained for similar weights in the case of the standard Laplacian in 1 10 1 . 
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